Testing the validity of the strong spin-orbit-coupling limit for octahedrally 
coordinated iridates in a model system Sr 3 CuIr06 
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The electronic structure of Sr3CuIr06, a model system for the 5d Ir ion in an octahedral environ- 
ment, is studied through a combination of resonant inelastic x-ray scattering (RIXS) and theoretical 
calculations. RIXS spectra at the Ir L3-edge reveal an Ir ti g manifold that is split into three lev- 
els, in contrast to the expectations of the strong spin-orbit-coupling limit. Effective Hamiltonian 
and ab inito quantum chemistry calculations find a strikingly large non-cubic crystal field splitting 
comparable to the spin-orbit coupling, which results in a strong mixing of the j e ff = | and j e ff = § 
states and modifies the isotropic wavefunctions on which many theoretical models are based. 



PACS numbers: 78.70.Dm,71.70.Ej,71.70.Ch,71.70.-d 

A new type of Mott physics has recently been discussed 
in the 5d transition metal oxides [TH3]- In contrast with 
the more familiar 3<i Mott insulators, for which new phys- 
ical phenomena originate from the large on-site Coulomb 
interaction, U, in 5d transition metal oxides the strong 
electron correlation has been argued to be driven by a 
large spin-orbit coupling (SOC). Iridates are an impor- 
tant example in this category and have attracted much 
attention recently [THUj- In particular, iridates in which 
the iridium atom sits in an octahedral oxygen cage have 
been classified as narrow-band Mott insulators. It is ar- 
gued that strong SOC splits the occupied t 2g orbitals 
into two sets of narrow bands corresponding to j e ff = 1 
and j e ff = | states. For the important case of the Ir 4+ , 
5d 5 ion, this puts one electron in the j e ft = \ state 
and even a small U, then gives rise to Mott insulating 
behavior. Further, because of the isotropic nature of 
the resulting wavefunctions [3], this has important conse- 
quences for the magnetic superexchange interactions and 
has led to predictions of an array of interesting physics, 
including quantum spin liquids [3 topological 
insulators [4] and superconductivity [14]. 

The strong SOC limit assumes local cubic symmetry 
with perfect IrOg octahedra. In reality, all of the pro- 
posed physical realizations of these models have non-ideal 
octahedra. Examples include the pyrochlores i?2 &2O7 
(R = Y, Sm, Eu and Lu), the hexagonal lattice sys- 
tems A 2 Ir03 (^4 = Li, Na), the hyperkagome lattice 
compound Na^^Os and the square lattice compounds 
(Ba,Sr) 2 Ir04. To date these have mostly been treated 
in the strong SOC limit. There are signs though, 
that this may not always be appropriate. For example 
Na2lr03 has been thought to be a realization of the Ki- 
taev model on a honeycomb lattice[5]- However, It was 



found recently to be antiferromagnetically ordered in a 
zigzag pattern p?CJl 12T] , contradicting theoretical predic- 
tions based on the strong SOC limit. This is suggestive 
of significant non-cubic crystal fields and indeed very 
recently theorists have begun to consider this question 
[221E3]. Thus, understanding the interplay of SOC and 
non-cubic crystal fields and experimentally determining 
their size are crucial requirements for any realistic model 
for the iridates. 




FIG. 1. (Color online) The Cu-Ir0 6 -Cu-Ir0 6 chain of 
Sr3CuIr06- The local coordinates for the IrOe octahedra and 
the CuOe prism are labeled with (x,y,z) and (x',y') respec- 
tively. 

In this work, we experimentally determine the spin- 
orbit coupling and crystal field splitting in a model iri- 
date, SraCuIrOglM]; which has only small distortions of 
the IrOe octahedra. We find that there is strong mix- 
ing of the j e ff = \ an d Jeff = § states and that, even 
in this case, it is not appropriate to work in the strong 
SOC limit. An important consequence of this mixing is 
that the wavefunctions are no longer isotropic and the 
magnetic superexchange interactions are modified away 
from the conventional picture. These results suggest that 
the conventional jeff = \ Mott insulator picture of many 
of the currently interesting octahedrally coordinated iri- 
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dates needs to be revisited in light of these effects. 

Before presenting the results in detail, we first note 
that there are several unique features of this study of 
SraCuIrOg that make it ideal for probing these effects. 
First, the oxygen octahedra are well separated from each 
other and therefore one does not need to consider hy- 
bridization effects between them [2 5], Second, the octa- 
hedra themselves only have relatively small distortions - 
the Ir-0 bond lengths vary by less than 2% (full details 
of the structure are provided below), so that non-cubic 
crystal fields are expected to be small. Third, we use the 
technique of resonant inelastic x-ray scattering to probe 
the energy and momentum dependence of the excitation 
spectrum. This technique has the advantage of prob- 
ing only the orbital and magnetic excitations associated 
with the Ir sites, making interpretation of the spectra 
straightforward. This, combined with the fact that the 
iridium-derived electronic states are well localized, allows 
us to compare to the theoretical models with substantial 
precision. We observe six well-defined electronic transi- 
tions. By comparing with an effective model Hamiltonian 
and with ab initio theoretical calculations, we are able 
to derive a complete description of the electronic degrees 
of freedom. 

SraCuIrOg is a quasi-one-dimensional material with a 
monoclinic structure [HI [25]. Individual IrOg octahedra 
are linked by spin- 1/2 Cu 2+ ions along one direction, 
forming a chain structure (fig. 1). Intersite hopping 
is largely suppressed due to this special structure: the 
d x 2_ y 2 orbital at the Cu site is essentially orthogonal 
to the t 2g orbitals on the Ir site, resulting in a negligible 
intersite hopping integral. Even though the extended na- 
ture of the 5d orbitals tends to reduce correlation effects, 
this special arrangement puts Sr 3 CuIr06 in the strongly 
localized regime. 
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FIG. 2. (Color online) Incident x-ray energy dependence of 
the RIXS spectra from SrsCuIrOe at the Ir L3 edge. Three 
low energy features at 0.28, 0.58 and 0.81 eV are observed. 
These show similar resonant behavior, each peaking at 11.216 
keV. 

Single crystal SraCuIrOg was grown with self- flux tech- 



niques. In reciprocal space, the chain direction is parallel 
to Q«(1.98 1). In our measurements, Q points at £(1.98 
l)+(0 13 0) were selected along the chain direction. A 
single chain is shown in Fig. 1 with the local coordinates 
for a single IrOg octahedra marked. The O-Ir-0 bonds 
are straight and of similar length with a variation of less 
than 2%. The distortion of the IrOg octahedra is mainly 
rotations of the O-Ir-0 bonds. The in-plane bonds along 
the x and y directions rotate towards each other, reducing 
the angle from 90 to 80 degrees, and the apical oxygens 
rotate towards the Cu atoms by ~ 4 degrees. We shall see 
that even these small distortions are sufficient to break 
the degeneracy of the t2 S manifold and move the system 
out of the strong SOC limit. 

The energy and momentum dependence of the excita- 
tion spectrum was studied using resonant inelastic x-ray 
scattering [27| . The experiments were carried out at the 
Ir L3 edge for which dipole transitions excite and de- 
excite a 2pa core-electron to the 5d orbitals. Due to 
2 

the strong core-hole potential and localized nature of the 
Ir 5d electrons in Sr 3 CuIr06, the RIXS process is dom- 
inated by local dipole transitions which lead to intra- 
site d-d excitations, an important consideration for this 
work. The measurements were carried out at beamline 
9-ID, Advanced Photon Source, in a horizontal scatter- 
ing geometry. A Si(844) secondary monochromator and 
a R=2m Si(844) diced analyzer were utilized. The overall 
energy resolution of this setup was ~ 45meV (FWHM). 
All data were collected at 7 K. 

The incident X-ray energy dependence of the RIXS 
spectra near the Ir L3 edge is shown in Fig. 2. Three fea- 
tures are observed at 0.28, 0.58 and 0.81 eV in the low 
energy region. Importantly, these three features show 
the same resonant behavior as function of incident x-ray 
energy, all resonating around 11.216 keV. This indicates 
that they all originate from initial 2p —> 5d transitions 
into the same unoccupied states within the Ir t 2g mani- 
fold. The different energies arise from the fact that the 
5d — > 2p de-excitation transitions create holes in different 
orbitals, leaving the system in different excited states at 
the end of the respective RIXS processes. 

More details of these three excitations can be seen in 
Fig. 3. The peak at 0.28eV, though sharing the same 
resonant behavior, has much less spectral weight. This 
excitation is most likely from impurity states associated 
with defects. At higher energy loss above 2eV, d-d exci- 
tations into Ir e g orbitals and charge transfer excitations 
from oxygen 2p orbitals are also observed. Here we focus 
on the 0.58 and 0.81 eV excitations, which are d-d exci- 
tation within the iridium t 2g manifold. These two peaks, 
especially the lower energy one, are fairly sharp. Impor- 
tantly, they show no significant dispersion with momen- 
tum transfer (Fig. 3) . These observations are consistent 
with the narrow band picture and the very localized na- 
ture of the 5d electrons in SraCuIrOg- 

In the ideal case for which the Ir atom sits within a 
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FIG. 3. (Color online) RIXS spectra at different Q = £(1-98 
l) + (0 13 0) points along the chain direction. The data at £ = 
-0.3 are fit with multiple Gaussian peaks. The black solid line 
shows the total fitting result. Individual Gaussian peaks are 
shown with red dashed lines, which are centered at 0.28, 0.58, 
0.81, 3.14, 3.90, 4.72 and 6.20 eV, respectively. The vertical 
and horizontal bars indicate the excitations listed in Table I 
from the MRCI+SOC calculations. 



regular oxygen octahedral cage, the ti g manifold is split 
into j e ff = | and j e ff = § states by the SOC[2E]. In 
this effective 2 level system, only one d-d transition is ex- 
pected. Our observations clearly deviate from this ideal 
picture, and even without any modeling, already point 
to the importance of non-cubic crystal fields. 

In order to understand these results, we first turn to 
an effective Hamiltonian to describe the local Ir 5d elec- 
trons. This will provide us with an intuitive description 
of the electronic structure. We will later see that this is 
consistent with a more sophisticated ab initio quantum 
chemistry calculation for this structure. We start with 
a local atomic multiplet Hamiltonian in hole representa- 
tion, 

(1) 

n no- 
where A is the strength of the SOC and A describes the 
splitting of the d xy from d yz ^ zx due to a non-octahedral 
crystal field, n is the index of the Ir sites. This Hamilto- 
nian splits the energy levels within the ti g manifold into 
three doublets, 

£0,2 = + 6 T V9 + 26 + P),E X = £ (2) 

where 5 = ^ and Eq is the ground state energy. Our 
RIXS spectra determine the energy differences to be E\ — 
E = 0.58 eV and E 2 - E = 0.81 eV, which leads to 
A = 0.44 eV and A = 0.31 eV. That is, the non-cubic 
crystal field is of a similar magnitude to SOC. With the 
parameters for Hi oca i in hand, the wavefunctions of the 



three doublets can now be determined. We find, 

10°) = rk(°- 65c W,; + *<W,T ± ^CTj.,t) 

\4>i) = ^(idyzu ± d ^,i.) (3) 

|02> = ^o&dxyt^ ~ 0.65(^4.,^ ± d zxi ^)] 

The above effective Hamiltonian description is simple but 
informative. The three doublet picture naturally explains 
the data with a SOC strength A that is consistent with 
earlier calculations pQ. Importantly, however, it suggests 
a significant non-cubic crystal field, one that is compa- 
rable with the SOC - despite the small distortions of 
the octahedra in this material. This non-cubic crystal 
field strongly modifies the relevant wavefunctions. For 
example, \4>q) shows significant deviation from the ideal 
case: instead of the expected equal admixture of the 
three t 2g orbitals in the ground state, the contribution 
from the d xy orbitals is reduced by more than 30% in 
SraCuIrOg. This significant modification is illustrated in 
Fig. 4. Given that anisotropic magnetic interactions lie 
at the heart of much of the proposed novel phenomena 
for the iridates, such a modification of the wavefunctions 
would be crucial and warrants further study. 




FIG. 4. (Color online) Density profile of the t2 g hole. Left: 
j e ff = 1 state for perfect IrOe octahedra. Right: the modified 
wavefunction |</>o) (Eqn. 3) due to the non-cubic crystal field 
arising from distortions of the octahedra. 

To do so, we performed ab initio many-body calcula- 
tions based on wavefunction electronic-structure theory. 
Complete-active-space self-consistent-field (CASSCF) 
and multircfcrence configuration-interaction (MRCI) 
techniques from modern quantum chemistry [29 , 3Qj were 
employed, as implemented in the MOLPRO package [31"]. 
Multiorbital and multiplet physics, SOC's, and O 2p to 
Ir hd charge transfer effects are all treated on an equal 
footing, fully ab initio. The calculations were carried out 
on a cluster which contains a central IrOg octahedron, 
two nearest-neighbor CUO4 plaquettes, and the adjacent 
Sr atoms, properly embedded in a large array of point 
charges that reproduces the crystal Madelung field in the 
cluster region. 

The effect of non-cubic crystal fields is studied by 
first turning off the SOC. Without SOC, the embedded- 
cluster MRCI calculations give a splitting of 0.45 eV be- 
tween the d XZ:VZ and d xy levels. This value is consistent 
with the A parameter derived in the modcl-Hamiltonian 
analysis carried out above. 
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TABLE I. The percentage contributions of the different Ir 
5a configurations to the different excitations. The results 
are obtained from MRC1+SOC calculations. 



Config. 



Energy(eV) 



0.59 0.91 3.22-4.84 5.88-5.91 



d 2 d 2 d 1 

^xy^xz^yz 

d 2 d 1 d 2 

11 x y ^xz^yz 
jl i2 i2 
^xy^xz^yz 
,4 1 
l 2g e g 

P 2 
z 2g e g 



37.1 49.1 13.8 
49.0 47.8 2.4 
13.9 3.1 83.8 



>96.0 



>99.0 



With spin-orbit interactions switched on 32J, the 
MRCI calculations predict excitations within the t^g 
manifold at 0.59 and 0.91 eV. At higher energies, we 
find Ir i| g to i|g e g excitations between 3.2 and 4.9 eV 
and O 2p to Ir 5d charge transfer transitions between 3.5 
and 5.4 eV. Further, excitation energies of ~5.9 eV are 
predicted for the Ir i!g e g states, see Table I. 

As seen in Fig. 3, distinct spectral features are found 
by RIXS in all those energy windows. The good agree- 
ment between the calculated values and the experimental 
observations confirms the reliability of applying MRCI 
calculations in this context. In Table I, the contribution 
of the different orbital configurations to each state within 
the manifold is listed. Again, those weights signifi- 
cantly deviate from the values corresponding to an undis- 
tortcd octahedron, in agreement with the effective model. 
For example, both approaches predict that in the ground- 
state configuration the hole is preferentially distributed 
over the d xz . yz orbitals, due to a sizable splitting between 
the d X z,yz and d xy levels. The ground-state wavefunction 
\4>o) from the effective model displays a hole occupation 
ratio of 0.65 2 :1:1, i.e., 0.18:0.41:0.41. This compares very 
well to the percentages found in the MRCI+SOC calcu- 
lations (see Table I, first column). The same holds for 
the other two doublets. 

Our experimental and theoretical studies in SraCuIrOe 
thus establish the importance of the local environment in 
determining the electronic structure of the Ir ti g mani- 
fold. With only a moderate distortion of the Ir0 6 oc- 
tahedra away from ideal cubic symmetry, the induced 
non-cubic crystal field has a similar strength to the SOC. 
This drives the system away from the strong SOC limit 
(-£ >> 1) and significantly modifies the relevant wave- 
functions, which will have important consequences for 
the magnetic exchange interactions. 

We now discuss the relevance of this work to other iri- 
dates, where to date, most of the theoretical work has 
assumed the strong SOC limit and treated the relevant 
electrons as in the effective j — ^ states in which d xyi 
d xz and d yz are equally mixed. Taking each of the iri- 
date families in turn, we find for the pyrochlores i?2 112O7 
that while all the O-Ir-0 bonds are straight and of equal 
length, they rotate towards each other by 6~10° [3"3"] . 



In Na^^Os, the O-Ir and Ir-0 bonds deviate from 90° 
by up to 9°, and differ in length by up to 5% [9]. For 
Na2lr03, the structure remains controversial. While a 
dramatic distortion was reported earlier [B], recent results 
suggest it is smaller with the O-Ir-0 bonds bending away 
from straight by ^7° and the O-Ir and Ir-0 bonds deviat- 
ing from 90° by up to 10° [21]. In all cases, the distortions 
are similar in magnitude or larger than the distortions 
present in Sr 3 CuIr0 6 . Thus treating these systems as be- 
ing in the strong SOC limit is likely to be incorrect. With 
the electronic structure and magnetic interactions highly 
dependant on the relative contribution from the different 
orbitals, the non-octahedral crystal fields must be explic- 
itly considered in any realistic models. Finally, we look 
at Sr 2 Ir04. In this system the octahedral bond angles are 
90°, and the bond lengths differ by only 4% [34 . Thus 
this system may be the closest to the strong SOC limit, 
consistent with experimental observations [28]. However, 
the present study suggests that this may be the exception 
rather than the rule. 

In summary, we have used resonant inelastic x-ray scat- 
tering, a microscopic model Hamiltonian and ab initio 
quantum chemistry calculations to study the electronic 
excitations in Sr 3 CuIr06, where Cu and IrOg octahedra 
form a 1-D spin chain. A three-level structure of the Ir 
t2 g manifold is observed and ascribed to the lowering of 
the local octahedral symmetry. We show that this lower- 
ing of symmetry leads to a strong non-octahedral crystal 
field, which is comparable to the spin-orbit-coupling and 
modifies the electronic structure, the ground state wave- 
functions and the magnetic exchange interactions signif- 
icantly. It should therefore be taken into account when 
modeling the magnetic behavior of all non-ideal iridates. 
Important examples in the current literature include the 
pyrochlores i?2lf207, the honeycomb lattice A2&O3 and 
the hyperkagome A 4 Ir 3 8 . 
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